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ABSTRACT 

Recent work on Fermi acceleration at ultra-relativistic shock waves has demonstrated the 
need for strong amplification of the background magnetic field on very short scales. Am- 
plification of the magnetic field by several orders of magnitude has also been suggested by 
observations of gamma-ray bursts afterglows, both in downstream and upstream plasmas. 
This paper addresses this issue of magnetic field generation in a relativistic shock precursor 
through micro-instabilities. In a generic superluminal configuration, the level of magnetiza- 
tion of the upstream plasma turns out to be a crucial parameter, notably because the length 
scale of the shock precursor is limited by the Larmor rotation of the accelerated particles in 
the background magnetic field and by the speed of the shock wave. We discuss in detail and 
calculate the growth rates of the following beam plasma instabilities seeded by the accelerated 
and reflected particle populations: for an unmagnetized shock, the Weibel and filamentation 
instabilities, as well as the Cerenkov resonant instabilities with electrostatic modes; for a 
magnetized shock, the Weibel instability and the resonant Cerenkov instabilities with the lon- 
gitudinal electrostatic modes, as well as the Alfven, Whisler and extraordinary modes. All 
these instabilities are generated upstream, then they are transmitted downstream. The modes 
excited by Cerenkov resonant instabilities take on particular importance with respect to the 
magnetisation of the downstream medium since, being plasma eigenmodes, they have a longer 
lifetime than the Weibel modes. We discuss the main limitation of the wave growth associ- 
ated with the length of precursor and the magnetisation of the upstream medium for both 
oblique and parallel relativistic shock waves. We also characterize the proper conditions to 
obtain Fermi acceleration at ultra-relativistic shock waves: for superluminal shock waves, the 
Fermi process works for values of the magnetization parameter below some critical value, and 
there is an intrinsic limitation of the achievable cosmic ray energy depending on the ratio of 
the magnetization to its critical value. We recover results of most recent particle-in-cell sim- 
ulations and conclude with some applications to astrophysical cases of interest. In particular, 
Fermi acceleration in pulsar winds is found to be unlikely whereas its development appears to 
hinge on the level of upstream magnetization in the case of ultra-relativistic gamma-ray burst 
external shock waves. 
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1 INTRODUCTION 

Substantial progress has been accomplished in this last decade on 
our theoretical understanding of the acceleration of particles at rel- 
ativistic shocks, revealing in more than one place crucial differ- 
ences with Fermi acceleration at non-relativistic shock waves. For 
instance, Gallant & Achterberg (1999), Achterberg et al. (2001) 
have emphasized the strong anisotropy of the cosmic ray popula- 
tion propagating upstream, which is directly related to the fact that 
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the relativistic shock wave is always trailing right behind the accel- 
erated particles. These particles are confined into a beam of open- 
ing angle 9 < l/r s h (with T s h the Lorentz factor of the shock 
wave in the upstream frame) and are overtaken by the shock wave 
on a timescale rL/F s h, with tl the typical Larmor radius of these 
particles in the background magnetic field. One consequence of 
the above is to restrict the energy gain per up — > down — > up cycle, 
AE/E, to a factor of order unity. Early Monte Carlo numerical ex- 
periments nonetheless observed efficient Fermi acceleration, with 
a generic spectral index s = 2. 2 — 2. 3 in the ultra-relativistic limit 
(Bednarz & Ostrowski 1998, Achterberg et al. 2001, Lemoine & 
Pelletier 2003, Ellison & Double 2004), in agreement with semi- 
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analytical studies (Kirk et al. 2000) and analytical calculations 
(Keshet & Waxman 2005). This value of the spectral index is how- 
ever restricted to the assumption of isotropic turbulence both up- 
stream and downstream of the shock (Niemiec & Ostrowski 2004; 
Lemoine & Revenu 2006), whereas the shock crossing conditions 
imply a mostly perpendicular magnetic field downstream, which 
severely limits the possibility of downstream scattering. Further- 
more, it was later stressed by Niemiec & Ostrowski (2006) and 
Lemoine, Pelletier & Revenu (2006) that these early studies im- 
plicitly ignored the correlation between the upstream and down- 
stream particle trajectories during a cycle. In particular, the former 
numerical study demonstrated that Fermi acceleration became in- 
efficient if the proper shock crossing conditions were applied to 
the background magnetic field. This result was demonstrated an- 
alytically in the latter study, concluding that Fermi acceleration 
could only proceed if strong turbulence (SB 3> B) existed on a 
scale much smaller than the typical larmor radius. The addition of 
turbulence on large scales ^> tl does not help in this respect, 
as the particle then experiences a roughly coherent field on the 
short length scales that it probes during its cycle. Further studies 
by Niemiec, Ostrowski & Pohl (2006) have confirmed that Fermi 
acceleration proceeds if short scale turbulence is excited to high 
levels, either downstream or upstream. The detailed conditions un- 
der which Fermi acceleration can proceed have been discussed ana- 
lytically in Pelletier, Lemoine & Marcowith (2009); they are found 
to agree with the numerical results of Niemiec, Ostrowski & Pohl 

(2006) . 

Amplification of magnetic fields on short spatial scales thus 
appears to be an essential ingredient in Fermi processes at ultra- 
relativistic shock waves. Quite interestingly, strong amplification 
has been inferred from the synchrotron interpretation of gamma-ray 
burst afterglows, downstream at the level of SB/ B > 10 4 — 10 5 
(Waxman 1997; see Piran 2005 for a review), and upstream with 
SB IB > 10 2 - 10 3 (Li & Waxman 2006), assuming an upstream 
magnetic field typical of the interstellar medium. Understanding the 
mechanism by which the magnetic field gets amplified is crucial 
to our understanding to relativistic Fermi acceleration, since the 
nature of this short scale turbulence will eventually determine the 
nature of scattering, hence the spectral index and the acceleration 
timescale. 

Concerning the amplification of the downstream magnetic 
field, the Weibel two stream instability operating in the shock tran- 
sition layer has been considered as a prime suspect (Gruzinov & 
Waxman 1999, Medvedev & Loeb 1999, Wiersma & Achterberg 
2004, Achterberg & Wiersma 2007, Achterberg, Wiersma & Nor- 
man 2007; Lyubarsky & Eichler 2006). Several questions never- 
theless remain open. For instance, Hededal & Nishikawa (2005) 
and Spitkovsky (2005) have observed, by the means of numerical 
simulations that this instability gets quenched when the magnetiza- 
tion of the upstream field becomes sufficiently large. On analytical 
grounds, Wiersma & Achterberg (2004), Achterberg & Wiersma 

(2007) , and Lyubarsky & Eichler (2006) have argued that it satu- 
rates at a level too low to explain the gamma-ray burst afterglow. 
The long term evolution of the generated turbulence also remains 
an open question, although Medvedev et al. (2005) claim to see the 
merging of current filaments into larger filaments through dedicated 
numerical experiments. 

Regarding upstream instabilities, the relativistic generaliza- 
tion of the non-resonant Bell instability has been investigated by 
Milosavljevic & Nakar (2006) and Reville, Kirk & Duffy (2006) in 
the case of parallel shock waves. However, ultra-relativistic shock 
waves are generically superluminal, with an essentially transverse 



magnetic field in the shock front. For this latter case, Pelletier, 
Lemoine & Marcowith (2009) have shown that the equivalent of 
the Bell non-resonant instability excites magnetosonic compressive 
modes and saturates at a moderate level SB/B ~ 1 in the frame 
of the linear theory. 

In recent years, particle-in-cell (PIC) simulations have become 
a key tool in the investigation of these various issues. Such simu- 
lations go (by construction) beyond the test particle approximation 
and may therefore probe the wave - particle relationship, which is 
central to all of the above issues. Of course, such benefice comes at 
the price of numerical limitations of the simulations, both in terms 
of dimensionality and of dynamic range, which in turns impact on 
the mass ratios accessible to the computation. Nonetheless, early 
PIC simulations have been able to simulate the interpenetration of 
relativistic flows and to study the development of two stream in- 
stabilities at early times, see e.g. Silva et al. (2003), Frederiksen et 
al. (2004), Hededal et al. (2004), Dieckmann (2005), Dieckmann, 
Drury & Shukla (2006), Dieckmann, Shukla & Drury (2006), 
Nishikawa et al. (2006), Nishikawa et al. (2007) and Frederiksen & 
Dieckmann (2008) for unmagnetized colliding plasma shells, and 
Nishikawa et al. (2003), Dieckmann, Eliasson & Shukla (2004a, 
b), Nishikawa et al. (2005) and Hededal & Nishikawa (2005) for 
studies of the magnetized case. The formation of the shock itself 
has been observed for both electron-positron and electron-proton 
plasmas thanks to recent simulations that were able to carry the in- 
tegration on to longer timescales, see e.g. Spitkovsky (2005), Kato 
(2007), Chang et al. (2008), Dieckmann, Shukla & Drury (2008), 
Spitkovsky (2008a, b), Keshet et al. (2009). All of the above studies 
use different techniques for the numerical integration, and varying 
parameters (dimensions, composition, mass ratios, density ratios 
of the colliding plasmas and relative Lorentz factors) in order to 
examine different aspects of the instabilities to various degrees of 
accuracy and over different timescales. 

Several of these studies have reported hints for particle ac- 
celeration through non Fermi processes (Dieckmann, Eliasson 
& Shukla 2004b; Frederiksen et al. 2004; Hededal et al. 2004; 
Hededal & Nishikawa 2005; Nishikawa et al. 2005; Dieckmann, 
Shukla & Drury 2006, 2008). Concrete evidence for Fermi acceler- 
ation, i.e. particles bouncing back and forth across the shock wave 
has come with the recent simulations of Spitkovsky (2008b), and 
was studied in more details for both magnetized and unmagnetized 
shock waves in Sironi & Spitkovsky (2009). In particular, this lat- 
ter study has demonstrated the inefficiency of Fermi acceleration 
at high upstream magnetization in the superluminal case, along 
with the absence of amplification of the magnetic field (thus in full 
agreement with the calculations of Lemoine, Pelletier & Revenu 
2006). This result is particularly interesting, because it suggests that 
the magnetization of the upstream plasma, in limiting the length 
of the precursor, may hamper the growth of small scale magnetic 
fields, and therefore inhibit Fermi cycles. Finally, the long term 
simulations of Keshet et al. (2009) have also observed a steady de- 
velopment of turbulence upstream of the shock wave, suggesting 
that as time proceeds, particles are accelerated to higher and higher 
energies and may thus stream further ahead of the shock wave. We 
will discuss this issue as well at the end of the present work. 

The main objective of this paper is to undertake a systematic 
study of micro-instabilities in the upstream medium of a relativistic 
shock wave. We should emphasize that we assume the shock struc- 
ture to exist and we concentrate our study on the shock transition 
region where the incoming upstream plasma collides with the shock 
reflected and shock accelerated ions that are moving towards up- 
stream infinity. Therefore, care should be taken when confronting 
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the present results to the above numerical simulations which repro- 
duce the collision of two neutral plasma flows in order to study the 
development of instabilities that eventually lead to the formation of 
the shock (through the thermalization of the electron and ion popu- 
lations). The physical set-up that we have in mind matches best that 
obtained in the simulations of shock formation and particle accel- 
eration described in Spitkovsky (2008b), Keshet et al. (2009) and 
Sironi & Spitkovsky (2009), or that simulated in Dieckmann, Elias- 
son and Shukla (2004a, b) and Frederiksen & Dieckmann (2008), 
or that studied in Medvedev & Zakutnyaya (2008). Our approach 
also rests on the following observation, namely that in the ultra- 
relativistic limit, the accelerated (or the reflected) particle popula- 
tion essentially behaves as an unmagnetized cold beam ofLorentz 
factor ~ r4. 

In the present paper, we assume the beam to be carrying a 
weak current and in so doing, we neglect electromagnetic current 
instabilities. We will nevertheless include in our summary of insta- 
bilities the relativistic generalization of the Bell current instability 
(Bell 2004), since it has been studied in detail in several recent stud- 
ies (Milosavljevic & Nakar 2006; Reville, Kirk & Duffy 2006). The 
instability triggered by the cosmic-ray current in the case of oblique 
shock waves has also been discussed in the relativistic regime in 
Pelletier, Lemoine & Marcowith (2009). Note also that in the case 
of pair plasmas, electromagnetic current instabilities do not take 
place as the beam remains neutral. 

The layout of the present paper is as follows. We examine the 
instabilities triggered by this beam, considering in turn the cases of 
an unmagnetized upstream plasma (Section 3) and that of a mag- 
netized plasma (Section 4). In Section 5, we discuss the intermedi- 
ate limit and construct a phase diagram indicating which instability 
prevails as a function of shock Lorentz factor and magnetization 
level. We then discuss the possibility of Fermi acceleration in the 
generated turbulence and apply these results to the case of gamma- 
ray bursts shock waves and pulsar winds. We will recover the trend 
announced above, namely that a magnetized upstream medium in- 
hibits the growth of the magnetic field hence particle acceleration. 
In Section 2, we first discuss the general structure of a collisionless 
shock, in the case of a electron-proton plasma with a quasi perpen- 
dicular mean field, borrowing from analyses in the non-relativistic 
limit. 



in the front frame). The electron temperature is expected to grow to 
a value comparable to, but likely different from that of protons, 



which reaches T„ 



(r sh - 1) 



The electrostatic barrier 



thus allows the reflection of a significant part of the incoming pro- 
tons since e$ ~ (r s h — 1) m p c . Although it reflects a fraction 
of protons, it favors the transmission of electrons that would other- 
wise be reflected by the magnetic barrier. The reflection of a frac- 
tion of the protons ensures the matter flux preservation against the 
mass density increase downstream. However because the magnetic 
field is almost transverse, an intense electric field E — fi s hB ener- 
gizes these reflected protons such that they eventually cross the bar- 
rier. Interactions between the different streams of protons are then 
expected to generate a turbulent heating of the proton population, 
which takes place mostly in the so-called "foot" region. This foot 
region extends from the barrier upstream over a length scale (in the 
shock front frame, as indicated by the | sh subscript) £ F | sh = r L | sh , 
where r L | s h denotes the Larmor radius of the reflected protons. 

Entropy production in the shock transition region comes from 
two independent anomalous (caused by collisionless effects) heat- 
ing processes for electrons and ions. The three ion beams in the foot 
(incoming, reflected in the foot and accelerated) interact through 
the "modified two stream instability", which seemingly constitutes 
the main thermalisation process of the ion population. A careful 
description of these anomalous heating processes certainly requires 
an appropriate kinetic description. For the time being, we note that 
the growth of the ion temperature develops on a length scale iy. 
The temperature of the electrons rather grows on a very short scale 
scale £r <C Iy which defines the "ramp" of the shock. In non- 
relativistic shocks, electrons reach a temperature larger than ions; 
however we do not know yet whether this is the case in relativis- 
tic shocks. These electrons also experience heating in the convec- 
tion electric field. Moreover, due to the strong gradient of magnetic 
field, an intense transverse electric current is concentrated, inducing 
anomalous heat transfer through the ramp. Probably an anomalous 
diffusion of electron temperature occurs that smoothes out the tem- 
perature profile; however it has not been identified in relativistic 
shocks. Electron heating is described by Ohm's law in the direction 
of the convection electric field (in the x x B direction, taken to be 
z): 



(1) 



2 GENERAL CONSIDERATIONS 

2.1 On the configuration of a relativistic collisionless shock 
wave 

A collisionless shock is built with the reflection of a fraction of 
incoming particles at some barrier, generally of electrostatic or 
magnetic nature. Let us sketch the general picture, borrowing from 
model of non-relativistic collisionless electro-ion shocks (see e.g. 
Treumann & Jaroschek 2008a, b for recent reviews). In an electron- 
proton- plasma carrying an oblique magnetic field, one expects a 
barrier of both electrostatic and magnetic nature to rise. Because 
the magnetic field is frozen in most part of the plasma, its trans- 
verse component is amplified by the velocity decrease. This in it- 
self forms a magnetic barrier which can reflect back a fraction of 
the incoming protons. Similarly, the increase of electron density 
together with the approach of the electron population towards sta- 
tistical equilibrium is concomittant with the rise of an electrostatic 
potential such that e$ ~ T e log(n/n u | sh ) (n is the local density in 
the front frame, and n u | sh the upstream incoming density viewed 



with fi x < in the shock front frame, E = /3 sh B u , B u denoting 
the background magnetic field at infinity. The magnetic field pro- 
file can be obtained by prescribing a velocity profile going from 
— fi s hC ~ — c to ~ — c/3 over a distance much larger than £r. The 
profile displays a ramp at scale £r followed by an overshoot before 
reaching the asymptotic value 3B U . The above result indicates that 
the relevant scale for £n is the relativistic resistive length: 



47T 



feff 



(2) 



This is a very short scale not larger than the electron inertial length 
5 e = c/uj vc even when the anomalous resistivity is so strong 
that the effective collision frequency v e g is of order uo pc . This 
scale thus represents the growth scale of three major quantities, 
namely, the potential, the magnetic field and the electron tem- 
perature. It is of interest to point out that this scale always re- 
mains much smaller than the foot scale. Indeed, even if S a is es- 
timated with ultra-relativistic electrons of relativistic mass r s h77i p , 
i.e. 5 e — [r s hm p c 2 /(47rn e | sh e 2 )] 1/2 , it remains smaller than the 
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foot length, since 



B \ B h 



47rn e | u r2 h r 



1/2 



< 1 



(3) 



using the value of £p| s h for particles with typical energy F s h m p c 2 
in the shock front. The last inequality in the above equation is a 
natural requirement for a strong shock. The downstream flow re- 
sults from the mixing of the flow of first crossing ions (adiabatically 
slowed down) with the flow of transmitted ions after reflection. All 
the ingredients of a shock are then realized. 

In the case of an electron-positron plasma, when a magnetic 
field is considered, no electrostatic barrier rises, only the magnetic 
barrier appears. However, if the mean magnetic field is negligible, 
a barrier can rise only through the excitation of waves, as demon- 
strated by the PIC simulations discussed above. 

The structure is thus described by two scales £r and If and 
three small parameters: £ cr , the fraction of thermal energy density 
behind the shock converted into cosmic ray energy, gb the ratio 
of magnetic energy density over the incoming energy density and 

1/Tah. 

2.2 Particle motion 

As mentioned above, there are three particle populations in the foot: 
the cold incoming particles, the reflected protons, and the acceler- 
ated particle population which has undergone at least one up^ 
down — > up cycle. This latter population arrives upstream with a 
typical Lorentz factor Tb ~ Fg h , with a typical relative spread of 
order unity. The second population of reflected protons also car- 
ries an energy ~ T 2 h m p c 2 , since these particles have performed a 
Fermi-like cycle, albeit in the front rather than downstream. There- 
fore one can treat these two populations as a single beam. From the 
point of view of the instabilities, one can approximate this beam 
as cold, with momentum distribution oc 5 (p x — F 2 h m p cj 5 (p±). 
Indeed, the instabilities are governed by the beam velocity, the 
dispersion of which remains very small, being of order A /3b ~ 
— (2/Fb)AFb/Fb. In order to verify this, one writes the suscepti- 
bility of the beam, assuming as above that it is unmagnetized on the 
scale of the instabilities (Melrose 1986): 



b 

X/ij 



Aire 



5u + 



/b(p) 

7 

kiC/3j + kjcfii (k 



k/3c 



+ 



2 )m 



,(4) 



with p = (3"/m p c and /b(p) the distribution function of the beam. 
Since the velocity distribution of the beam is essentially delta like, 
one may then indeed approximate the above beam susceptibility 
with that of a cold beam; the difference amounts to a redefinition 
of the beam plasma frequency by a factor of order unity. 

Another crucial length scale in our study is the length scale of 
the precursor. As discussed above, this length scale £p| s h = r L | sh 
in the front shock in the case of a magnetized shock wave. In the 
upstream frame, this can be rewritten as: 



rL\u 

1 sh 



cJciFsh sin Ob 



(Bu ^ 0) 



(5) 



We assume that the field is almost perpendicular in the front 
frame, but in the upstream comoving frame we consider its obliq- 
uity (angle 6b with respect to the shock normal), assuming that 
sin 6b > l/r s h. The particular case of a parallel shock wave for 
which r s h sin 6b < 1 is discussed in Section [B~2l there it will be 



shown that a fraction (1 — r a h#s) of the particles that return up- 
stream may actually propagate to upstream infinity in the limit of a 
fully coherent upstream magnetic field, while Eq. {5) remains cor- 
rect for the rest of the accelerated particle population. The size of 
the precursor for the particles that escape away is eventually given 
by the level of turbulence ahead of the shock wave. 

In the case of an unmagnetized shock wave, the size of the pre- 
cursor is determined by the length traveled by the reflected protons 
in the self-generated short scale turbulence. Neglecting for simplic- 
ity the influence of the short scale upstream electric fields (we will 
see in Section [6] that this does not affect the following result), this 
length scale can be written (Milosavljevic & Nakar 2006; Pelletier, 
Lemoine & Marcowith 2009): 



r L|u 

F 4 / 



(6) 



where £ c represents the typical scale of short scale magnetic fluc- 
tuations. Whether one or the other formula applies depends on sev- 
eral possible situations and outcomes: if the shock is magnetized 
and one considers the first generation of cosmic rays, one should 
use Eq. lO; if the shock is magnetized and one assumes that a sta- 
tionary state has developed with strong self-generated turbulence, 
one should use Eq. (|6); obviously, if the development of the turbu- 
lence cannot take place, one should rather use Eq. l|5}; finally, for an 
unmagnetized shock, Eq. ([6} applies. In the following, we discuss 
the turbulence growth rate for these different cases. 

There seems to be a consensus according to which magnetic 
fluctuations have to be tremendously amplified through the gener- 
ation of cosmic rays upstream in order for Fermi acceleration to 
proceed. A fraction £ cr of the incoming energy is converted into 
cosmic rays and a fraction of this cosmic rays energy is converted 
into electromagnetic fluctuations, which add up to a fraction £ e m of 
the incoming energy. This process is expected to develop such that 
the generation of cosmic rays allows the generation of electromag- 
netic waves that in turn, through more intense scattering, allows 
further cosmic ray acceleration and so on until some saturation oc- 
curs. We write the quantities f cr and £ em as: 

f = — P = ^ cm n\ 

r sh n u m p c T^ h n u m p c ! 

with £ om < f cr . We approximate the beam pressure with that of 
the cosmic rays, i.e. P CI ~ r s hn b | s h m p c 2 for the first generation 
of accelerated particles, as expressed in the shock front frame. The 
electromagnetic energy density is written I7 cm in the same frame, 
as usual. 

Unless otherwise noted, our discussion takes place in the up- 
stream rest frame in what follows. 



3 UPSTREAM INSTABILITIES IN THE ABSENCE OF A 
MEAN MAGNETIC FIELD 

When the ambient magnetic field can be neglected or is absent, the 
reflected particles and the fraction of particles that participate to the 
first Fermi cycle constitute a relativistic cold beam that pervades the 
ambient plasma and trigger three major micro-instabilities. One is 
the two stream electrostatic instability, which amplifies the elec- 
trostatic Langmuir field through a Cerenkov resonant interaction 
uj — k ■ Vb = 0, with k || E || Vb. Another is the Weibel instabil- 
ity, with k || Vb -L E and its analog filamentation instability, with 
k _L v b || E (Bret, Firpo & Deutsch 2004, 2005a, 2005b; see also 
Bret 2009 for a recent compilation). These two instabilities are non- 
resonant and mostly electromagnetic with a low phase velocity so 
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that the magnetic component of the wave is dominant. It is thus par- 
ticularly relevant for developing particle scattering. Finally, these 
authors have also discovered an oblique resonance which grows 
faster than the above two. It is mostly longitudinal (see further be- 
low) but k is neither perpendicular nor parallel to the beam. These 
growth rates are easily recovered as follows. 

For a cold beam, Eq. l[4j gives the following susceptibility: 



b _ 

Xij 



"pb 



Su + 



kiC/3bj + kjc[3 hl (k 2 c 2 - u 2 )(3 h if3bj 



k-/3 b c 



(u-k-(3 b cY 



(8) 

The beam propagates with velocity /3bC = (l — l/r^) 1 2 x; the 
relativistic beam plasma frequency (in the upstream frame) is given 
by: 



I^pb 



/ 47rn b | u e ; 

V r b m p 



1/2 



(9) 



recalling Fb — Fg h . One can solve the dispersion relation, includ- 
ing the beam response, to first order in x° since its contribution is 
of order: 



Wpb\ _ me_ 

Up c J 



Ccr < 1 



(10) 



hence a growth rate: 



V3 

2 4 /3 



1/3 



(16) 



One should note that the Cerenkov resonance can only take place 
with plasma modes with phase velocity smaller than c (refraction 
index kc/ui(k) > 1), hence transverse modes are excluded in this 
respect. 

The oblique mode, with k z 7^ and a resonance as above 
yields a growth rate that is larger by a factor T^ 3 than the two 
stream rate given in Eq. J16t for k z — (Bret, Firpo & Deutsch 
2004, 2005a, b). This can be understood as follows. The instability 
arises from the xx component of the beam susceptibility tensor, 
which dominates over the other components at the resonance [see 
Eq. GO], and which reads: 



b 

/Cxx 



u 2 h u 2 /r 2 + pjklc 2 

u 2 (co - p h k x cf 



(17) 



This component is suppressed by 1/T b when k z — 0, which ex- 
plains the factor appearing in the r.h.s. of Eq. U6b . For fc z 7^ 
however, the algebra is more cumbersome. Nevertheless, proceed- 
ing as above, with the resonance condition Eq. dl4t . one obtains in 
the limit 5 < 1 and ^ 1: 



5.0.7 Weibel/filamentation instability 

Consider now a mode with k y = 0, but k x 7^ 0, k z 7^ 0. The 
dispersion relation, including the beam response can be written as 
follows, to first order in %ij ■ 

2 2 , 2 2 b 2\ 

XyyU ) 



2 7 2 2 
■Wp-K C 



( 2 2 , 2 2 b W 2 2 , 2 2 . b \ 

[u -u p - k z c + x xx J [u -u p - k x c + \ zz j 



( 



k x k z c 2 + x xz u 2 



= . 



(11) 



with ui = + u 2 e . In the limit k x — > 0, one recovers the 
filamentation (Weibel like) instability by developing the above dis- 
persion relation to first order in x b , with: 



2 

-^ P b 



Up + k 2 c 2 



(12) 



It saturates at a growth rate J(uwc 



Upb in the limit kc 3> u p . 



3.0.2 Cerenkov resonance with oblique electrostatic modes 

In the other limit k z — » 0, one can simplify the dispersion relation 
for electrostatic modes down to: 



2 2 . b 2 

w - Up + x xx u 



. 



(13) 



Then, the two stream instability resonance condition between the 
Langmuir modes and the beam reads: 



2k 2 



(18) 



In the limit k z ^> k x 
oblique mode: 



Up/c, one recovers the growth rate of the 



T(u) 



V3 / 2 \V3 
^J73 [UpbUp) 



(19) 



This mode obviously grows faster than the previous two. 

Obviously, the mode is quasi-longitudinal, since resonance 
takes place with the electrostatic modes. However it also comprises 
a small electromagnetic component, | B y \ / \ E z \ ~ 2 1<5[, as can be 
seen by solving for the eigenmode, using the full dispersion relation 
including the beam contribution. 



4 INSTABILITIES IN THE PRESENCE OF A MEAN 
FIELD 

As before, we look for an instability of the upstream plasma waves, 
triggered by the beam of accelerated (and shock reflected) parti- 
cles. At non-relativistic shocks, one usually considers an interac- 
tion at the Larmor resonance. However this cannot be relevant in 
the ultra-relativistic case, because the interaction must develop on 
a distance scale < If which is itself much shorter than the Lar- 
mor radius. The particular case of a relativistic parallel shock wave 
will be briefly discussed thereafter. Note finally that for the fre- 
quently valid condition f3 A^sh sin < 1, the precursor has a 
length much larger than the minimum scale for MHD description 
(^mhd/^p| u = /3a F s h sin 9b), which justifies the resonance be- 
tween the beam and the MHD modes. 



u = u p (1 + 5) = Pbk x c(l + 5) 



(14) 



with by assumption |5| <C 1. After insertion into Eq. dl3t . this 
yields: 



2 

^pb 

2I>! 



(15) 



4.1 Oblique magnetic field 

In order to excite fast waves of frequency higher than the Lar- 
mor frequency, we consider again the Cerenkov resonance be- 
tween the non-magnetized beam and the magnetized plasma waves: 
u — k • Vb = 0. Let us recall that for a ultra-relativistic beam, 



6 M. Lemoine and G. Pelletier 



the velocity distribution is strongly peaked at ~ c, even if the 
dispersion in Lorentz factor of the beam is significant. We also dis- 
cuss the possibility of generating the magnetic field through a (non- 
resonant) Weibel (filamentation) instability with k x = 0. 



4.1.1 Weibel - filamentation instability 

This instability taking place in the shock transition layer between 
the unshocked plasma and the shocked plasma has been dis- 
cussed in detail in the waterbag approximation for an unmagnetized 
plasma (Medvedev & Loeb 1999; Wiersma & Achterberg 2004; 
Lyubarsky & Eichler 2006; Achterberg & Wiersma 2007, Achter- 
berg, Wiersma & Norman 2007). As we now argue, the Weibel in- 
stability can also proceed in the regime of unmagnetized proton - 
magnetized plasma electrons at smaller frequencies, corresponding 
to the range u c i < w < w cc (see also Achterberg & Wiersma 
2007). Again, we should stress that we consider a pure ion beam 
(reflected and accelerated particles), whereas most above studies 
consider two neutral interpenetrating plasmas. 

To simplify the algebra, we write down the dispersion relation 
in a frame in which the (x, z) plane has been rotated in such a way 
as to align B with the third axis, denoted zs; y remains the second 
axis ys . To simplify further a cumbersome algebra, we consider a 
wavenumber k || ys, perpendicular to both the beam motion and 
the magnetic field. The plasma di-electric tensor is written in this 
B frame as: 



£l - 1J 
— i£2 





IE2 
El 




(20) 



with the following usual definitions (for uj c i <ti u> <§C u; cc ): 



£i 



2 2 
Upj Ul pe 

LO 2 CO 2 



£2 



LO 2 



and 7] = kc/ui. One needs to rotate the beam susceptibility tensor 
to this B frame. The quantity of interest will turn out to be the 3 — 3 
component Xz 

X° , the dispersion relation indeed has the solution: 



= cos 2 8bXxx + sin 2 9bXzz- To first order in 



2 . 2 

r) + cos 



9s xL + sin 2 I 



'BXz 



. 



(22) 



Given the dependence of Xxx on ^> this is a quartic equation which 
admits the solution leading to the Weibel (filamentation) instability: 



2 2 a kc 



(23) 



As in the unmagnetized case, it saturates at a growth rate ~ 
o; p b cos 9b (up to the angular dependence on B). Note that in the 
limit cos 9b — * 0, this instability does not disappear. In order to see 
this, one has to consider the other branch of the dispersion relation, 
for cos 9b = 0, k = k z z: 



(24) 



One of the roots corresponds to the Whistler mode and the other to 
the Weibel unstable mode with to 2 ~ — u> 2 b . 

The above thus shows that fast waves can be excited by the 
relativistic stream in the intermediate range between MHD and 
electron dynamics, i.e. with unmagnetized plasma ions but magne- 
tized electrons. The typical length scale of these waves for which 
maximal growth occurs is obviously the electron inertial scale 
S c = c/ujp as before. 



4.1.2 Cerenkov resonance with longitudinal modes 

The previous discussion of the Cerenkov instability with electro- 
static modes can be generalized to the magnetized plasma limit by 
considering those modes with k || B, which do not feel the mag- 
netic field (see Lyubarsky 2002 for a discussion of this instability 
in the case of pulsar magnetospheres). Rewriting the above plasma 
di-electric tensor for a wavenumber parallel to the magnetic field, 
it is straightforward to see that the dispersion relation admits the 
longitudinal branch given by: 

£|| + cos 2 9 B Xxx + sin 2 9 B x\z = . (25) 

In order to avoid confusion, it may be useful to stress that the 
previous oblique denomination refers to the angle between the 
wavenumber and the beam direction, while the present term lon- 
gitudinal here refers to the parallel nature of k and B. At the 
Cerenkov resonance uo = lu p (1+5) = f3bk x c(l+S), with \5\ <§C 1, 
one has \Xxx I ^> \ Xzz I > therefore one can obtain the following ap- 
proximate solution for the growth rate: 

, "I 1/3 



V3_ 

2 4 /3 



2 2 a 

W pb ljJp COS 6b 



1 



Recalling that k x = k cos 9b = lo p /c, and that V\ 
can neglect in all generality the first term in the parenthesis, so that 



(26) 
sh, one 



(wpb^p 



.1/3 



4.1.3 Resonant instability with Alfven modes 

Turning now to resonant instabilities with Alfven waves, we con- 
sider a wavector in the (x, z) plane. The resonance condition for 
Alfven modes reads: k x ~ /3a k cos 9k, where 9^ represents the 
angle between the wavenumber and the magnetic field direction. 
Since (3 a < 1, this implies k x <€. k, therefore the wavenumber is 
mostly aligned along z and 9^ ~ tv/2 — 9b- 

The plasma dielectric tensor now reads (we omitted negligible 
contributions in sin 2 9k): 



A 



£i — rf cos 2 9k 

— 1£2 



l£2 



J) COS t>k Sin t>k 





ij\B = I -«£2 £1-7) 

\ 7j cos 9k sin 9k en— r/ 2 sin 2 9 k 



with (ll) < cj ci ): 



£i 



1 n w p 

-2- , £2 ~ , £|| ~ 1 

Pa u 2 



, (27) 



(28) 



The beam susceptibility can be approximated accurately by ne- 
glecting all components in front of Xx X , which dominates at the 
resonance, as explained above. The relevant components then are: 



b 

Xx B x B 
b 

Xx B z B 



■ 2 a b b 2 a b 

sin 9bXxx , Xz B z B - cos 9 B Xxx , 
Xz B x B — sin9s cos9bX xx • 



(29) 



The dispersion relation then takes the form: 

-j- - k 2 c cos 2 9 k ) (w 2 + k 2 c sin 2 9 k ) 
Pa J 

+k 4 c 4 sin 2 9 k cos 2 9 k - cu 4 A xx x xx = , (30) 

where A xx ~ - sin 2 9 B ^l/to 2 in the limit fc<5 e < 1. Writ- 
ing down the resonance condition uj — f^Ak cos 9kC (1 + 5) = 
(3bk x c (1 + 5), with \S\ <§C 1 as before, one obtains the growth 
rate: 



- ^3 (^lbPAkccos9 k ) 1/ ' i 



(3D 
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where we approximated k z ~ k; recall furthermore that cos 8k — 
sin 8b ■ This instability disappears in the limit of a parallel shock 
wave as one can no longer satisfy the Cerenkov resonance condi- 
tion. 

One should stress that the above perturbative treatment re- 
mains valid as long as the condition \5\ <C 1, which amounts to 
£cr <C 0x at maximum wave growth rate (fee = oj c i). Therefore 
Alfven growth is limited to strongly magnetized shock waves only. 

In the continuity of right Alfven waves (the left modes being 
absorbed at the ion-cyclotron resonance), there are Whistler waves 
for quasi parallel propagation (with respect to the mean field), that 
are electromagnetic waves with a dominant magnetic component. 
For quasi perpendicular propagation, there are the ionic extraordi- 
nary modes, which have frequencies between the ion-cyclotron fre- 
quency and the low-hybrid frequency (obtained for large refraction 
index) and which are mostly electrostatic with a weaker electro- 
magnetic component. For scattering purpose, the whistler waves 
are the most interesting in this intermediate range; they are ac- 
tually excited in the foot of non-relativistic collisionless shocks 
in space plasmas. But for pre-heating purposes, the extraordinary 
ionic modes are more interesting (they are actually used for addi- 
tional heating in tokamaks). Let us now discuss these in turn. 



4.1.4 Resonant instability with Whistler waves 

We proceed as before, using the plasma di-electric tensor Eq. Ml\ 
in the range u> C i < w < lo cc with the components given in 
Eq. J2U . The Whistler branch of the dispersion relation reads, to 
first order in the beam response x° approximated by Eq. \29\ : 



2 2 

ti — r\ cos I 



h + Xxx i 



n ) -4 = . (32) 



When the beam response is absent, one recovers the dispersion re- 
lation for oblique Whistler waves: 



2 w cc , 4 4 2 a 

w wh . ~ — — k c cos 8 k . 



(33) 



Introducing the resonance lo = LOwh. (1 + 6) = Pbk x c (1 + 5), 
with \8\ <C 1, we obtain the growth rate: 



T( \ V^3 / 2 \l/3 



(34) 



In the latter equation, we again approximated k z ~ k, since the 
resonance condition implies k x <C k (therefore cos 8 k ~ sin^s). 
The instability disappears in the limit of a parallel shock wave as 
well, because the resonance condition cannot be satisfied. Maxi- 
mum growth occurs here as well for k ~ c/u pc ~ c/ui v , i.e. at 
the electron inertial scale S c , however the excitation range extends 
to the proton inertial scale 5i where it matches with the Alfven wave 
instability. 

As before, the perturbative treatment remains valid as long as 
\S\ <C 1, which amounts to £ CI <C (m p /m e ) 2 f)\. This condition 
is more easily satisfied that the corresponding one for amplification 
of Alfven waves. It will be discussed in more detail in Section|5] 

4.1.5 Resonant instability with extraordinary modes 

At MHD scales, the extraordinary ionic modes (that propagate with 
wave vectors almost perpendicular to the magnetic field) assimi- 
late to magneto-sonic modes. These modes has been shown to be 
unstable when there is a net electric charge carried by the cosmic 
rays (Pelletier, Lemoine & Marcowith 2009). The obtained growth 



rates are increasing with wave numbers indicating an instability that 
reaches its maximum growth at scales shorter than the MHD range. 
Let us therefore discuss how this instability extends to sub-MHD 
scales. 

Let us first discuss the ionic (lower hybrid) branch, lo < wn, 
with wih = t/u>^uJ^. In the following, we assume for simplicity 
lo cc <C u>pe, i-e- a weakly magnetized plasma. In the B frame, in 
which B is along zb and the beam propagates in the (x, z) plane, 
take k || y_g, with a small component k XB , i.e. in the (x, z) plane 
but perpendicular to B. The dispersion relation to zeroth order in 
X b reads: 



2 e 
V = - 



with (since lo < u>ih <§C lo cc ): 



4 



4 



C-l 



4 



(35) 



(36) 



hence 






c 2 f 2 

e l e 2 

ei 






e l e 2 




2 

W pc 






ui 2 h - LO 2 



(u c i «W « L0 lh ) 



(37) 



At lo <C uJci, this gives the fast magnetosonic branch with loh 
P^kc, while at lo c i <C lo <§C com, loh 
We define: 



T>(k,u) = 



ei 



4 



wmfcc/ ^/k 2 c 2 + oo 2 c . 
2 ■ (38) 



so that: 



2 _d_ 

1 duo 2 



V(k,uj) 



2 

71 ~H 



2 



~ rj (u <C w c i) , 

(U) c i <g LO <g LO\h) 



(39) 



Including the beam response, the dispersion relation becomes: 



2 
«2 



t\T) + (ei 



2 



i sin 2 ObxIx = . (40) 



We neglect the term r/ 2 B <C t] 2 in front of t\ ~ 1/0A ( at u 
LOci)- At the resonance lo — uih(1 + 5), with loh the solution of 
T>(k,iLiu) = 0, one finds: 



1 Wpb 
2 



2 sin 2 8 B 



Lo 2 -^-V(k,to) 

Old- 1 



(41) 



The growth rate for Cerenkov resonance with the lower hybrid ex- 
traordinary mode thus reads: 



T(w L x) 
Z(colx) 



2 4/3 

V3_ 

2 4/3 



Lj 2 b sin 2 8 B -r^PAkc 



1/3 



k 2 

2 ■ 2 a K y 

tjpbsm 8b 



uimUJ pe kc 



(lo <g lo c 

1/3 



k2 (k 2 c 2 + lo 2 c ) 3/2 _ 

(LO c i <C LO <ti L0 lh ) . 



(42) 



In the limit of magnetosonic modes, lo <§C w c j, one recovers 
the same growth rate as for Alfven waves; note that /3a kc <S 
o; c i implies k <C w p ;/c. At smaller scales, one finds that the 
growth rate reaches its maximum at k ~ lo pc /c with T(wlx) 

/ 2 

cient heating of the protons in the foot. 



sin 2 dBLOih) 1 ^ 3 ■ We can expect this instability to provide effi- 
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c 2 c 2 

£l — t'2 



Turning to the electronic (upper hybrid) modes, around uj 



with uj x 

(ujl + UJ, 



2 \l/2 



Wpe(^ 2 - U>\ 
LOpc — CJ ce /2, OJ z 



(43) 



uj pc + lu cc /2 and uj u h 



The dispersion relation takes the same form 



T>(k, uj) = 0, but now: 

a / 2 2 2 

O _,, , 9 / UJ UJ UJ 

—V(k,uj) ~ r/ 2 — j + — j - — + 1 

OljJ 1 \uj a — uj x UJ — u>i UJ — UJ h 

(44) 

The growth rate can be written in the same algebraic form as <4U . It 
vanishes in both limits uj — *• uj x and uj — > uj z , while for uj ~ uj pc , 
giving r\ ~ 1, one obtains: 



TV ~i ~ ^ ^2 . 2 fl ^ 2 c ^ 



(45) 



It vanishes in the limit uj cc /uj pc — > 0, in which limit the electronic 
extraordinary branch actually disappears. 

Being electrostatic in nature, these waves participate mostly 
to the heating process in the shock foot or precursor. However their 
scattering efficiency is comparable to the magnetic perturbations as 
will be seen further on. 



4.2 The particular case of a parallel magnetic field 

When the magnetic field is almost parallel, i.e. 9b < 1/T s h, the 
relativistic Bell non-resonant instability (Bell 2004, 2005) can de- 
velop (e.g. Milosavljevic & Nakar 2006; Reville, Kirk & Duffy 
2006). This instability is triggered by the charge current carried 
by the cosmic rays in the precursor, which induces a return current 
in the plasma, thereby destabilizing non-resonant waves of wave- 
length shorter than the typical Larmor radius, the cosmic rays being 
unresponsive to the excitation of the waves. The growth rate of this 
instability in the upstream frame is (Reville, Kirk & Duffy 2006): 



1 (WBell) 



A>n b i u 



(46) 



and growth is maximal at the scale k c ~ T(uiBen) /(Pac). 

One can then verify that, under quite general assumptions, this 
growth rate is larger than the growth rate of the Weibel instability, 
since the ratio of these two is given by: 



J(^Bcll) 
T (uJWe. ) 



r 2 .1/2 
shScr 



(47) 



One must emphasize however that the Bell instability is quenched 
when the growth rate exceeds the ion cyclotron frequency see 
Couch, Milosavljevic & Nakar (2008), Riquelme & Spitkovsky 
(2009), Ohira et al. (2009). This limitation will be made clear in 
Section|5T2l 



5 LIMITATIONS OF THE INSTABILITIES 

Using the growth rates derived previously, we can now delimit the 
conditions under which the various instabilities become effective, 
and which one dominates. We then discuss the limit between un- 
magnetized and magnetized shock waves, from the point of view 
of these upstream instabilities. 



5.1 Superluminal shock waves 

In this Section, we discuss the generic case of relativistic superlu- 
minal shock waves, taking sin 2 9 b ~ 1. Unless otherwise noted, 
we assume an e — p plasma; we will discuss how the results are 
modified in the limit of a pair plasma at the end of this discus- 
sion. The more particular case of relativistic parallel shock waves 
is treated further below. 

We start by introducing the two parameters X and Y defined 
as follows: 



X 
Y 



1 sh , 

m p 

f>2 

r 4 -°u|u 
1 sh" 



47rnb| u m p c 2 



r 2 f -l 
1 sh^u^cr 



(48) 



The upstream magnetization parameter a u also corresponds to 
the Alfven velocity squared of the upstream plasma (in units of 
c 2 ). If the field is fully perpendicular, the shock crossing con- 
ditions imply Bdld.J- — 5 u | u ,xr s h\/8, and for the enthalpy 
/i<j|d — (8/3)r 2 h /i u | u (for a cold upstream plasma, see Bland- 
ford & McKee 1976), so that cr d ~ 3<t u sin 2 6>b. If the mag- 
netic field is mostly parallel, meaning sin#g ^ l/r s h, then 
a A ~ (3/8)r- h V u . 

Let us first compare the growth rates of the instabilities ob- 
tained in the magnetized case; the unmagnetized case (in particular 
the oblique mode) will be discussed thereafter. We carry out this 
comparison at the wavenumber where the growth rates reach their 
maximum, namely k ~ uj pc /c. The ratio of the Weibel to Whistler 
instability growth rates is given by: 



1 (uwc.) 



1/6 



(49) 



hence the Weibel instability will dominate over the Whistler 
Cerenkov resonant instability whenever Y <C X. The fastest 
mode however corresponds to the Cerenkov resonance with 
the longitudinal modes along the magnetic field, since the ra- 
tio of the growth rates of this mode to the Weibel mode is 
(m p /m e ) 1/ ' 6 ^CT 1 ^ 6 , which is always greater than one. 

Since the Cerenkov resonant instabilities for the Whistler and 
Alfven waves scale in a similar way with the eigenfrequencies of 
the resonant plasma modes, it is straightforward to see that Whistler 
waves will always grow faster than the Alfven waves. 

Concerning the extraordinary modes, one finds that 
2T(a;wh.)/I(u'Lx) ~ {rrip/m^) 1 ^' on the ionic (lower hy- 
brid) branch, while I(wwh.)/I(wux) ~ {uJ-pc/uJcc) 1 ^ 3 on the 
electronic (upper hybrid) branch. Therefore the growth of these 
modes is always sub-dominant with respect to that of Whistler 
and Weibel modes. Since the growth rates of the Alfven and 
extraordinary modes are always smaller than that of the Whistler 
modes, we discard the former in the following. 

Additional constraints can be obtained as follows. First of all, 
the above derivation of the instabilities has assumed the beam to 
be unmagnetized, i.e. that the growth time be much shorter than 
the Larmor time of the beam particles. This condition is always 
easily satisfied, since it reads: Y <C rf h for the Weibel instabil- 
ity, Y <C rf h ^ cr 1/ ^ 3 (m e /m p ) _1 ' /3 for the longitudinal mode and 
Y ^ rf h m p /m e for the Whistler Cerenkov resonant mode. One 
can explicit the dependence of Y on the shock parameters in or- 
der to verify this; for the Weibel instability, the condition amounts 
to £„■ ^> T~L <r n , which is indeed easily verified at large Lorentz 
factors. 

Concerning the Cerenkov resonant modes, one must also re- 
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quire that |<5| < 1 in order for the perturbative treatment to be 
apply. In the case of longitudinal modes, this is automatically sat- 
isfied since u p b < u p . However, for Whistler modes of smaller 
eigenfrequency, this implies a non-trivial constraint ui p ^ < aj-yvh. 
which can be translated into Y ^> X 2 for wwh. = <^cc- 

Further bounds can be obtained by requiring that the back- 
ground protons are non-magnetized in the case of the Weibel in- 
stability, which requires X (lo) 3> w<». This condition is however 
superseded by the requirement that the growth can occur on the 
precursor length scale, since £f/c ~ (r s hWci) _1 [see Eq. l|5}]. 
At this stage, it is important to point out a fundamental difference 
between the Cerenkov resonant instabilities and the Weibel / Ala- 
mentation instabilities. The former have, by definition of the res- 
onance, a phase velocity along the shock normal which, to zeroth 
order in \5\ exceeds the shock velocity, while the latter have van- 
ishing phase velocity along x. Therefore the timescale available for 
the growth of these non-resonant waves is the crossing time of the 
precursor: they are sourced at a typical distance £f away from the 
shock, then advected downstream on this timescale. Regarding the 
resonant modes, their phase velocity along x is 0^^ = /3b(l + <5j?), 
with 8r — TZ(S). Since Sr < for the resonant modes, one must 
consider three possible cases: (i) /3$ lX < Ash, in which case the 
mode is advected away on a timescale £f /c as for the non-resonant 
modes; (ii) j3$ iX > /3 s h, in which case the mode propagates for- 
ward, but exits the precursor (where it is sourced) on a similar 
timescale; and (iii) j3$ }X ~ /3 s h, in which case the mode can be 
excited on a timescale ~ c~ 1 £f/( Ah — A> } x ) and where the diver- 
gence corresponds to the situation of a mode surfing on the shock 
precursor. However condition (i) appears to be the most likely, as 
least in the ultra-relativistic limit, for it amounts to 2T 2 h | <5r j 3> 1 . 
Indeed, all resonant instabilities have a growth rate ~ (lo^uj) 1 ^ 3 
where to is the eigenfrequency of the resonant mode (an excep- 
tion is the upper hybrid mode for which the growth rate is smaller 
by (ajcc/ajpc) 2 / 3 , in which case the following condition is even 
stronger), therefore the condition 2r 2 h |<5fj,j ^> 1 can be rewritten 
as w/up C < 7(r sh /10) 3 (£cr/0.1) 1/2 , which is generically sat- 
isfied. This means that the phase velocity of the resonant modes, 
when corrected by the effect of the beam becomes smaller than 
the shock front velocity, so that these modes are advected on a 
timescale ~ £f/c and transmitted downstream, after all. For the 
purpose of magnetic field amplification downstream and particle 
acceleration, this is certainly noteworthy, as such true plasma eigen- 
modes (Whistler, Alfven, extraordinary and electrostatic longitudi- 
nal or oblique modes) can be expected to have a longer lifetime 
than the Weibel modes. 

The modes thus grow on the precursor crossing timescale if 
T (uj) If /c 2> 1, which can be recast as 7 < 1 for the Weibel 
instability, Y < £™c?(m e /m p )~ 1/3 and XY < 1 for the 
Cerenkov resonant Whistler mode. Henceforth, we use the param- 
eter G = (Lu ph /cj p )- 2/:i =£ c -r 1/3 (m e /m p )- 1/3 > 1. 

In short, we find that the various instabilities discussed here 
are more likely quenched by advection rather than by saturation. 
In Section loTTl we provide several concrete estimates for cases of 
astrophysical interest and it will be found that this limit is indeed 
quite stringent. 

Finally, one must also require that the growth rate of the 
Cerenkov resonant instabilities does not exceed the proper eigen- 
frequency of the mode. For the Whistler modes, as discussed at the 
end of Section l4.1.4l this implies Y > X 2 . 

In Section 3, we have also examined the growth rates in 
the absence of a mean magnetic field, and concluded that the 
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Figure 1. Instability diagram for superluminal relativistic shock waves (as- 
suming sin 2 8g ~ 1): in abscissa, X = T s ^m e /m v , in ordinates Y = 

Fj h B 2 / (47rn b | u m p c 2 ). The parameter G = C c " r 1/3 (m e /m P )" 1/3 > 
1. The axes are plotted in log-log on arbitrary scale. The main result is 
summarized by the thick solid line, which indicates the maximum value 
of Y(X) which allows electromagnetic waves to grow. The other lines 
indicate the boundaries of the regions of growth of the various instabili- 
ties, as indicated. The hierarchy of growth rates, from largest to smallest 
is as follows: oblique and longitudinal, then Whistler and/or Weibel. The 
long dashed line separates the regions in which the growth of Whistler or 
Weibel modes is faster: for values of Y(X) larger than the long dashed 
line, Whistler modes grow faster. The growth rates of the oblique mode 
and the longitudinal mode are comparable. Unlike the longitudinal mode, 
the oblique instability is limited by the assumption of unmagnetization (see 
main text), but at the same time, it applies to a larger wavenumber phase 
space. The regions for Alfven and extraordinary modes are not indicated 
(see main text). 

oblique mode of Bret, Firpo & Deutsch (2004, 2005a, b) was by 
far the fastest. This instability is very similar to the Cerenkov res- 
onance with the longitudinal modes propagating along the mag- 
netic field and indeed the growth rates only differ by sin 2 / 3 9b, see 
Eqs. <19M26t. The difference lies in the degree of magnetization 
of the ambient plasma: while the oblique mode is limited to the un- 
magnetized limit, the longitudinal mode does not suffer from such 
constraint; the oblique mode however covers a larger fraction of the 
wavenumber phase space than the longitudinal mode. 

With respect to the oblique mode instability, the shock can be 
described as unmagnetized as long as the background electrons and 
protons remain unmagnetized on the timescale of the instability; of 
course, one must also require that the instability has time to grow on 
the length scale of the precursor. Note that the latter condition also 
implies that the beam can be considered as unmagnetized over the 
instability growth timescale, which is another necessary condition. 
For the oblique modes, those conditions amount to: 

X(woW.) > ^ & Y < GX 2 , (50) 
2T(woM.) > c/If Y < G, (51) 

with G = £CT 1 ^ 3 (m e /m p ) -1/ ' 3 > 1 as above. Provided the 
above two conditions are satisfied, the oblique mode dominates 
over the Weibel and Whistler Cerenkov instability growth rates, 
just as the longitudinal mode. The Cerenkov resonant instabil- 
ity with Whistler waves dominates over the oblique modes when 
X < G _1/3 and GX 2 < Y < l/X. The Cerenkov reso- 
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nant instability with Whistler waves dominates over the longitudi- 
nal modes when X < G^andG < Y < 1/X.ForX < G" 1 
and Y > X' 1 , or for G" 1 < X and Y > G neither of the 
above instabilities can grow. For reference, X < G _1 corresponds 
to r 8h < 150fc/ 3 . The above regions can be summarized in the 
X — Y plane as in Fig. [TJ which delimit the domains in which 
the various instabilities can grow, and which of these instabilities 
dominates in each case. 

From the above discussion, the case of a pair shock is eas- 
ily obtained by taking m p /m e — > 1, by restricting oneself to the 
study of the oblique and Weibel modes, and by considering only the 
right part of Fig. [TJ with X > 1, since X = r s h for a pair shock. 
One sees that, irrespectively of F s h, the oblique and longitudinal 
modes can grow if Y < £ C r and the Weibel mode grows on the 
precursor timescale if Y < 1, 

In this respect, it is instructive to compare the present re- 
sults with the latest simulations of Sironi & Spitkovsky (2009). 
These authors find that the growth of instabilities is quenched when 
the magnetization <r u > 0.03 for a perpendicular (or oblique) 
pair shock with T s h — 20. This corresponds to X ~ 20 and 
Y ~ lO^ 1 (<7 u /0.03). Our results indicate that indeed, at this 
high level of magnetisation, both Weibel and oblique/longitudinal 
instabilities are quenched by advection. Note that these simulations 
do not exclude that the instabilities are quenched even at lower 
magnetisations. Our calculations thus bring to light the following 
point of interest. One should not infer from the simulations of 
Sironi & Spitkovsky (2009) that superluminal shock waves can- 
not lead to magnetic field amplification. This conclusion entirely 
depends on the level of magnetization. It would therefore be inter- 
esting to extend the PIC simulations down to weakly magnetized 
shocks with a u ~ 3 x 10 _3 £ 2 / 3 (r s h/20)~ 2 in order to probe the 
limit at which the oblique mode can grow. 



5.2 Parallel shock waves 



In the ultra- relativistic limit, parallel shock waves are non-generic; 
however, they may lead more easily to particle acceleration than su- 
perluminal shock waves (since the argument discussed in Lemoine, 
Pelletier & Revenu 2006 no longer applies) and consequently pro- 
vide interesting observational signatures. One can extend the above 
discussion to the case of parallel shock waves as follows. 

First of all, the main limitation of the instabilities, that is 
due to the precursor crossing timescale disappears in the limit 
r s h sin 9b — * as a fraction 

1 — /3 sh cos 9 B — tJ— sin Ob „ 
P = -=-2* ~ (1 - r sh 9 B ) 2 (52) 

1 — Psh 

of the particles can propagate to upstream infinity (at least in the 
limit of a fully coherent magnetic field). This can be seen as fol- 
lows. Particles cross the shock wave back toward downstream once 
their angle cosine with the shock normal becomes smaller than /3 s h- 
However, when F s h sin 6b < 1, there exists a cone Cb around the 
magnetic field direction, of opening angle 9 B — acos(/3 s h), that 
never intersects the cone C defined around the shock normal with 
opening angle acos(/3 s h) . Because the pitch angle of the particles 
with respect to the magnetic field direction is conserved, particles 
that enter toward upstream in this cone Cb never recross the shock 
toward downstream (up to the influence of the turbulence). The 
fraction of particles that enter in this cone is approximately given 
by the ratio of the solid angles, i.e. the factor (1 — T s ^9b) 2 quoted 
before. Depending on the value of r s h#B, this fraction can be sub- 




m e /m, 



Figure 2. Same as Fig.[T]for the case of a parallel relativistic shock wave. 
The regions of growth are drawn at a fixed value of r s h sin 9b (here taken 
to be 0.3). The axes are plotted in log-log on arbitrary scale. The longitudi- 
nal mode can grow in all parameter space (in the idealized limit of a fully 
coherent upstream magnetic field) due to the divergence of the precursor 
length (see main text). 



stantial and these particles can excite plasma waves up to large dis- 
tances from the shock. Of course, the actual precursor length re- 
mains finite as a result of the influence of large and short scale tur- 
bulence. In the following, we will simply discard these advection 
constraints in order to avoid introducing new parameters. 

We also choose to discuss the limitations as a function of 
X and Y, but at a fixed value of sin6*s < 1. The fact that 
X oc r s h and that T s h sin 9b is fixed modifies slightly the limita- 
tions derived previously. The main limitation for the oblique mode 
is the non-magnetization condition, X(w bi. ) 3> w cc which can be 
rewritten Y <C GX 2 as before. For the Weibel mode, the con- 
dition of non-magnetization of the protons, I(ww e . ) 3> uJd now 
reads Y < (m p /m e ) 2 X 2 . 

The Cerenkov resonant mode with longitudinal waves does 
not suffer any constraint in this parallel configuration, but the 
growth rate now becomes significantly smaller than that of the 
oblique mode (when the latter applies), by a factor sin 2 ' 3 9b- 

Concerning the Whistler modes, the condition \S\ < 1 now 
leads to Y 2> (m e /m p ) _2 (r s h sin#.g) _2 X 4 . One recovers the 
condition expressed in the case of superluminal shock waves for 
sin 9% ~ 1, as one should. At a fixed value of r s h sin 0b however, 
the condition appears slightly different. Note that in the particular 
case of parallel shock waves, there is another non-trivial condition 
for these Whistler modes, which is related to the fact that the eigen- 
frequency u>wh. — w C e sin 9b (at maximal growth) should exceed 
oj c i, in order for the Whistler branch to apply. This translates into 
X < r sh sin0B < 1. 

Finally, the Bell non-resonant instability requires the back- 
ground protons to be magnetized, as noted in Section |4~2l which 
corresponds to Y 2> £cr (rn p /m e ) 6 X 6 . 

These various conditions are expressed in Fig. [2] which is the 
analog of Fig.QJfor parallel shock waves. Here as well, the limit of a 
pair shock can be obtained simply by taking the limit m e /m p — > 1 
and discarding the Whistler branch as well as the Bell instability 
(which requires a net current to exist upstream). 

As for the oblique shock wave, it is instructive to compare the 
present results to the simulations of Sironi & Spitkovsky (2009), 
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who find in particular that instabilities can be triggered in parallel 
shocks for a magnetisation ass = 0.1. Here ass corresponds to 
the definition of the magnetisation given in Sironi & Spitkovsky 
(2009), or to our defintion of downstream magnetisation up to a 
factor 3/4, the latter factor of 3/4 representing the difference be- 
tween enthalpy and energy density for a relativistic gas. For a par- 
allel shock wave, this thus corresponds to an upstream magneti- 
sation cr u = 0.05lt h , hence to Y ~ O.OSlt^" 1 ~ lO 4 ^" 1 
for X ~ 20. One may note that at this large level of magneti- 



sation, one has lu c . 



One can check immediately, using 



the above, that neither the oblique nor the Weibel instabilities can 
grow, as their respective non-magnetisation conditions are not sat- 
isfied. The growth of instabilities observed in the simulations of 
Sironi & Spitkovsky (2009) is thus likely due to the Cerenkov res- 
onance with the longitudinal modes that propagate along the mag- 
netic field, which does not suffer from these constraints. 

The above also allows to understand the abrupt transition as 
a function of shock obliquity: when sin 6 b — ► 1/Tsh, the growth 
of fluctuations is suddenly inhibited and so is Fermi acceleration. 
Indeed, as sin 6 b —* l/r s h, the fraction of particles that can es- 
cape to upstream infinity vanishes, hence the precursor length now 
rapidly decreases to the value given by Eq. (|5}. This prevents the 
growth of fluctuations at high magnetisation levels, as discussed 
above for superluminal shock waves, and consequently this pre- 
vents successful Fermi acceleration (Lemoine, Pelletier & Revenu 
2006). 



6 TRIGGERING FERMI ACCELERATION 

It is important to underline that Fig.[T]indicates whether instabilities 
triggered by the first generation of cosmic rays returning upstream 
have time to grow or not. If these instabilities cannot be triggered 
by the first generation, meaning if the shock wave characteristics 
are such that (X, Y) lie above the thick solid line of Fig.Q] then 
instabilities cannot be triggered, either upstream or downstream (at 
least in the frame of our approach). Consequently Fermi cycles will 
not develop, at least for superluminal shock waves, in accordance 
with the arguments of Lemoine, Pelletier & Revenu 2006, Pelletier, 
Lemoine & Marcowith 2009 and with the simulations of Niemiec, 
Ostrowski & Pohl 2006. 

If, however, the initial values of X and Y are such that in- 
stabilities can develop, Fig. Q] suggest that these instabilities will 
develop upstream and be transferred downstream. Fermi cycles 
may then develop provided the appropriate conditions discussed in 
Lemoine, Pelletier & Revenu (2006) and Pelletier, Lemoine & Mar- 
cowith (2009) are satisfied. These conditions have been discussed 
under the assumption of isotropic short scale magnetic turbulence, 
and we restrict ourselves to this assumption in the present work 
as well. It would certainly be interesting to generalize this discus- 
sion to more realistic turbulence configurations, as in Hededal et al. 
(2004), Dieckmann, Drury & Shukla (2006) for instance. However, 
this clearly becomes more model dependent in terms of turbulence 
configuration and for this reason, we postpone such a study to fu- 
ture work. 

Let us discuss first the case of upstream turbulence. When par- 
ticles are scattered off short scale £ c , but intense magnetic fluctua- 
tions, the scattering frequency of a relativistic particle of momen- 
tum p is 



Since the oblique mode dominates over the Whistler and Weibel 
waves over most of the parameter space, one cannot ignore the in- 
fluence of short scale electrostatic fields. These electrostatic waves 
lead to a second order Fermi process in the upstream medium, with 
a concomittant pitch angle scattering. Indeed, the particle scatters 
against random electric fields iEu along the shock normal (x di- 
rection), gaining momentum Apii — i&EiiAt, with At ~ cj" 1 
at each interaction, and similarly in the perpendicular direction. The 
initial pitch angle of the particle (with respect to the shock normal) 
6 <C 1 in the upstream frame, and the particle is overtaken by the 
shock wave whenever 6 > l/r s h (Achterberg et al. 2001). This 
pitch angle diffuses according to: 

<A0 2 > ^ <A P 2 > ^ ,El + 26^E\ 

-T c , (54) 



At 



p 2 At 



for a correlation time r c = £ c /c ~ ui~ c . Therefore we obtain a 
scattering rate similar to the previous one d53b in which the mag- 
netic field fluctuation is replaced by the electric field fluctuation: 



(55) 



This correspondence justifies that we treat the short scale elec- 
tric and magnetic fields on a similar footing and consider the to- 
tal electromagnetic energy content. A conversion of a fraction of 
the energy of the beam into magnetic or electrostatic fluctuations 
is expected with £ cm < £ cr , with typically f cr ~ 10 _1 and 
Com ~ 10~ 2 - 10" 1 (Spitkovsky 2008a). Scattering in the short 
scale electromagnetic turbulence will govern the scattering process 
if itleadsto {Ap 2 )/p 2 ~ l/r 2 h on a timescale r L | S /(r sh c), with 
r L | B the Larmor radius of first generation cosmic rays as measured 
upstream relatively to the background magnetic field (see the cor- 
responding discussion in Pelletier, Lemoine & Marcowith 2009). 
If this short scale turbulence governs the scattering process, then 
Fermi acceleration will operate. Assuming £ c = c/cu pc , this condi- 
tion amounts to: 

1/2 



fe, 

Using the fact that £ c 
bound on a u : 



>r sh (^ 



1/2 



(56) 



< £ cr , this constraint can be rewritten as a 



cr u < Ccr l ah 



(57) 



This limit is very stringent indeed; in terms of our above param- 
eters, it can rewritten as F < X£ cr /r a h. We will discuss the 
applicability of this inequality in concrete cases in the following 
sub-section. 

If this condition is not verified, the background unamplified 
magnetic field remains the main agent of particle scattering up- 
stream. In this case, Fermi acceleration cycles can develop only 
if short scale turbulence govern the scattering downstream of the 
shock wave. As discussed in Pelletier, Lemoine & Marcowith 
(2009), this requires: 



8B ld 

4|d < r L id < s — 4|d 

-Did 



(58) 



(53) 



where all quantities should be evaluated in the downstream rest 
frame, and r L | d refers to the Larmor radius of the accelerated par- 
ticles in this frame. This double inequality amounts to requiring 
fhatf c | d /c < r s < 7l,o, i.e. that the scattering time r s = v^ 1 be 
shorter than the Larmor time in the mean field tl,o in order to break 
the inhibition constraint of the mean field that tends to drag the par- 
ticles in the downstream flow. The scattering must also develop in 
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a special regime where the correlation time £ c \d/c is shorter than 
the Larmor time. Regarding £ c | d , two main spatial scales are to be 
envisaged: the previous upstream electron skin depth, if one as- 
sumes that the typical scale of transverse fluctuations is preserved 
through shock crossing, and the downstream electron skin depth, 
if reorganization takes place through shock crossing. Assuming a 
typical electron temperature ~ r s hm p c 2 behind the shock, and ac- 
counting for shock compression of the electron density, this latter 
scale can actually be written as c lo ~. (lo pi the upstream ion plasma 
frequency), a factor 43 larger than the previous one. One should 
also envisage the possibility that the turbulence spectrum evolves to 
larger scales with time (Medvedev et al. 2005; Lemoine & Revenu 
2006; Katz, Keshet & Waxman 2007) but we will not do so here. 
Let us consider the above two possibilities in turn. 

If ^ c |d = c/ojpc (upstream electron skin depth), then the first 
inequality in Eq. l |58t can be rewritten as £ cm < m p /m e and is 
therefore always satisfied. The second inequality amounts to <7d < 
(m e /m p )^ 2 m , hence Y < r sh X£ 2 m /£ cr . This latter inequality 
is much more stringent. If satisfied, it means that the downstream 
short scale turbulence governs the scattering process, in particular it 
allows the particle to escape its orbit around the shock compressed 
background magnetic field on a timescale smaller than the Larmor 
time in this field. This is a necessary condition for successful Fermi 
cycles. 

If ^c|d = CLU pi (equivalently, the downstream electron skin 
depth), then the first inequality in Eq. d58fc becomes £ cm < 1, 
which is always true. The second inequality reads Od < £ 2 m ( or 
Y < r 2 h £ 2 m /£ cr ). We will summarize the two above two possible 
cases for i c \& and parameterize the uncertainty on £ g m by writing 
the condition as: 

(7 d < (J» = K Com , (59) 

with £ c | d = nc/ujpi and m e /m p < k < 1. One should how- 
ever recall that the typical scale of electromagnetic fluctuations 
could evolve with the distance to the shock front, as envisaged 
in Medvedev et al. (2005), Lemoine & Revenu (2006) and Katz, 
Keshet & Waxman (2007). This amounts to making k be a growing 
function of the energy taking values larger than 1. The above re- 
sult clearly reveals the need for dedicated PIC simulations of shock 
wave at moderate magnetisation, with realistic proton to mass ratio 
and geometry in order to reduce this large uncertainty on k and de- 
termine the precise conditions under which Fermi acceleration can 
take place. 

To summarize this discussion, we obtain the following condi- 
tions for successful Fermi acceleration. If Eq. l |57t is satisfied [or, 
to be more accurate, Eq. l|56H, then Fermi acceleration will op- 
erate, because the short scale fluctuations produced upstream are 
sufficiently intense to govern the scattering. In this case, it is im- 
portant to stress that Eq. @, which defines the length of the pre- 
cursor, no longer applies. It should be replaced by Eq. {6), which 
is larger. Physically, the precursor widens, giving more time for the 
fluctuations to grow, thus reaching a higher efficiency in terms of 
£em/£cr- If Eq. l !57t is not satisfied, e.g. because the upstream mag- 
netization is not small enough, particles gyrate in the background 
magnetic field before experiencing the short scale turbulence. Then 
Fermi acceleration will operate if Eq. d59t is verified. Conse- 
quently, a sufficient condition for the development of Fermi cycles 
is £ em > (<7 cri t/fi;) 1/2 , or equivalently a, > a clit , where a clit 
is the maximum value of the upstream magnetization that allows 
turbulence to grow upstream and then be transferred downstream. 
As shown previously for a superluminal configuration, a CI i t = 
F~ h ^ CT m v jm R for an electron-proton plasma in the realistic case 



where T s < 150£cr so that the transition is governed by the ex- 
citation of whistler waves; or <r cr it = r,^ 2 ^ 2 / 3 (m p /m e ) 1//3 for 
an electron-positron plasma, with the development of the oblique 
two stream instability. The spectral index and the maximal en- 
ergy remain to be determined however. In this respect, we note 
that Eq. l !58t provides an upper bound for this maximal energy: 
e max ~ r s hm p c 2 (a»/(j u ) 1 ^ 2 in the front frame. 

The more likely development of the Fermi process is thus hy- 
brid, in the sense that it is of drift type upstream and of diffusive 
type downstream. As Fermi cycles develop, particles are acceler- 
ated beyond the energy r 2 h m p c 2 considered here for the first gen- 
eration. Although they are less numerous, they stream farther ahead 
of the shock and are therefore liable to induce stronger amplifica- 
tion. One can only speculate about these issues, since the spectral 
index depends strongly on the assumption made on the shape of the 
turbulence spectra, upstream as well as downstream. In particular, 
if the magnetic field amplified downstream through the Weibel in- 
stability decays on scales of order of tens or hundreds of electron 
inertial lengths <5 C , the particles will likely escape towards down- 
stream because of the lack of scattering agents, thereby cutting 
off the Fermi process prematurely. Nevertheless, assuming for the 
sake of discussion that Fermi cycles develop with a spectral index 
s ~ 2 — 3, the number density of cosmic rays streaming upstream 
scales as n b | u (> p,) oc (p„/po) 1-s , with po ~ r 2 h m p c 2 . The 
beam plasma frequency, which controls the growth rates of the 
instabilities, w p »(> p») oc (p,/po) _S//2 , whereas the precursor 
length ^p| u (> p-t) oc (p,/po). Since the growth rates of the res- 

2/3 

onant instabilities which develop upstream scale as uj pt , s < 3 
would guarantee that the growth factor of the instabilities triggered 
by these high energy particles exceeds that for the first generation. 
These findings seem in agreement with the numerical simulations 
of Keshet et al. (2009) and Sironi & Spitkovsky (2009) who observe 
wave growth farther from the shock from high energy particles, as 
time increases. 

6.1 Applications 

It is interesting to situate the relativistic shock waves of physical 
interest in the above diagram. Here we consider three proto-typical 
cases: a pulsar wind, a gamma- ray burst external shock waves ex- 
panding in the interstellar medium, and a gamma-ray burst external 
shock wave propagating along a density gradient in a Wolf-Rayet 
wind. We find the following: 

• Pulsar winds: with T ~ 10 6 and a u ~ 0.01, one finds 
(X, Y) ~ (500, lO 10 ^" 1 ); the level of magnetization is thus so 
high that no wave can grow, either upstream or downstream. Fermi 
acceleration should consequently be inhibited. 

• Gamma-ray burst external shock waves expanding in the inter- 
stellar medium: for Y ~ 300 and a u ~ 10" 9 (i.e. B ~ 3 /xG), 
one finds (X, Y) ~ (0.1, lO" 5 ^. 1 ). Wave growth should be ef- 
ficient both usptream and downstream. Concerning Fermi acceler- 
ation, Eq. J57t amounts to Y < ( cr m e /m p . It can thus be only 
marginally satisfied. However, Eq. l |59t is most likely satisfied, so 
that Fermi acceleration should develop, even in the early afterglow 
phase when T s h ~ 300. 

• Gamma-ray burst external shock waves propagating along a 
density gradient in a Wolf-Rayet wind: taking a surface magnetic 
field of 1000 G for a WRq Wolf-Rayet progenitor, the magnetiza- 
tion at distances of 10 17 cm is a u ~ 10 -4 (Crowther 2007). This 
gives (X, Y) ~ (0.1, Ccr 1 )- Growth may or may not occur in this 
case, depending on the precise values of F s h, a u and £ cr . In detail, 
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the condition for Weibel growth Y < 1 is likely not verified for 
the above fiducial values, but could be verified in less magnetized 
winds and at later stages of evolution, with a smaller value of T s h. 
The condition for growth of Whistler waves, Y < l/X, may be 
satisfied if £ cr > 0.1 and it is likely to be more easily verified at 
smaller values of r s h and a u . Finally, the (most stringent) condi- 
tion for growth of the oblique mode, Eq. l !50t . is likely not verified 
in the initial stages with r s h ~ 300 and the above fiducial value of 
(T u , but would be verified if a u were smaller. 

However, Eq. J57b cannot be satisfied in this case, meaning that 
the orbit of the particle upstream is governed by the wind magnetic 
field, not by the amplified short scale component. Regarding the 
bound Eq. d59t , it can be satisfied, depending on the values of the 
wind magnetisation and most particularly on the value of k. The 
possibility of Fermi acceleration thus remains open in this case. 
More work is necessary to understand the properties of downstream 
turbulence in order to determine whether particle can eventually be 
accelerated. 



6.2 Further considerations 

It is important to emphasize that we do not understand yet the 
structure of a relativistic shock front in detail. In the previous 
section we have assumed that the shock front is structured like 
a non-relativistic front and just extended the non-relativistic re- 
sults. Since MHD compressive instability and extraordinary ionic 
modes can be excited, we cannot exclude that the foot be full 
of relativistically hot protons and electrons of similar tempera- 
ture •ynipC 2 with 1 <C 7 ^ T s . In that case the plasma re- 
sponse would be different, because the intermediate whistler range 
(and also extraordinary range) would disappear so that the plasma 
would behave like a relativistic pair plasma. Then, the relevant 
instabilities are the Weibel and oblique modes (in the unmagne- 
tized approximation). The length of the precursor and the Weibel 
growth rate remain unchanged, hence the domain of growth of 
the Weibel instability also remains unchanged. The growth rate 
of the oblique mode is however reduced because the background 
plasma frequency is smaller by a ratio (7m p /m e ) 1 ' 2 . Therefore 
the condition of growth on the advection timescale now reads 
Y < ^ 1/3 r^ 3 X- 1/3 (7m p /m e )" 1/3 . The ratio of the growth 
rates of the oblique mode to the Weibel mode can be written as 
(7?cr) -1,/6 , hence the Weibel instability becomes the dominant 
mode if 7 2> . 

In the downstream plasma, the magnetic fluctuations gener- 
ated by the Weibel instability are expected to disappear rapidly be- 
cause they do not correspond to plasma modes. Whistler and other 
resonant eigenmodes (when they are excited) are however transmit- 
ted and although they are not excited downstream, their damping is 
weak. When Fermi cycles develop, they create "inverted" distribu- 
tion downstream, that should produce a maser effect. 

Tangled magnetic field carried by the upstream flow are very 
compressed downstream and thus opposite polarization field lines 
come close together. This should produce magnetic reconnections 
in an unusual regime where protons and electrons have a similar 
relativistic mass of order F s h m p c 2 . Such a regime of reconnection 
deserves a specific investigation with appropriate numerical 
simulations. Despite magnetic dissipation, reconnections would 
probably create a chaotic flow that favors diffusion of particles 
from downstream to upstream. 



7 CONCLUSIONS 

In this work, we have carried out a detailed study of the micro- 
instabilities at play in the precursor of a ultra-relativistic shock 
wave. The main limitation for the growth of these waves is related 
to the length of precursor, which is itself related to the level of 
magnetisation in the upstream plasma (where magnetisation refers 
to the background field, not the shock generated short scale fields). 
Nevertheless, we have found electronic and ionic instabilities that 
grow sufficiently fast in the precursor of a relativistic shock. The 
fastest growing instabilities are due to the Cerenkov resonance be- 
tween the beam of accelerated (and shock reflected protons) and 
the upstream plasma Whistler waves and electrostatic modes. The 
Weibel instability, which is non-resonant by essence, is also ex- 
cited, but its growth is generally superseded by that of the previ- 
ous modes. The strongest amplification occurs on very short spatial 
scales ~ S c , the electron skin depth in the upstream plasma. Our 
results are summarized in Fig.[T]for the generic case of relativistic 
superluminal shock waves, which delimits the domains in which 
electromagnetic modes are excited in terms of shock Lorentz fac- 
tor and upstream magnetisation and defines the critical value <T cr M 
of the magnetization below which Fermi process can operate. Fig- 
ure|2]presents the corresponding limitations for the case of parallel 
shock waves; in this case, the growth of instabilities is made much 
easier by the divergence of the precursor length for a fraction of the 
particles returning upstream. Our results explain some features of 
recent PIC simulations of relativistic pair shocks of various geome- 
tries and magnetisation levels. 

We have discussed the conditions under which Fermi accel- 
eration can proceed superluminal shock waves once a significant 
fraction of the cosmic ray energy has been dumped into these short 
scale electromagnetic fluctuations. Fermi acceleration can operate 
if the upstream magnetisation (<r u ) or downstream magnetisation 
(<7d) are low enough for the shock generated turbulence to gov- 
ern the scattering of particles. This is the second condition that 
states the required level of electromagnetic energy density versus 
the magnetization. This requires either a u <C Com( m e/m-p)r~ h 2 
(for upstream scattering), which is however difficult to fulfill or 
o"d <SC «£em (f° r downstream scattering, which is easily ful- 
filled with a low level of turbulence); £ cm indicates the fraction 
of incoming energy transferred into electromagnetic fluctuations, 
with fern ~ 10~ 2 — 10 _1 generally indicated by PIC simula- 
tions, and k is a fudge factor that encaptures our ignorance of the 
transfer of electromagnetic modes excited upstream through the 
shock, m e /m p < « < 1 (and it may be even larger depend- 
ing on the particle energy if the scale of the electromagnetic flu- 
ctuations evolves with the distance to the shock). We emphasize 
the need for PIC simulations with realistic geometry, realistic pro- 
ton to electron mass ratios and moderate magnetisation (of order of 
the above) in order to lift this uncertainty on k and to determine 
the precise conditions under which Fermi acceleration can take 
place. This limitation also places a strict upper bound on the max- 
imum energy that is achievable through the Fermi process, namely 
e max ~ r sh m p c 2 ((T cr it/o"u) 1/2 , with a cri t the maximum magneti- 
sation that allows waves to grow, and <r u the upstream background 
magnetisation (see discussion after Eq. l |59l l. Beyond this intrinsic 
limit, the scattering time indeed becomes longer than the Larmor 
time in the mean field downstream, so that the particle is adverted 
downstream by the mean field and Fermi cycles end. 

We have also applied our calculations to several cases of astro- 
physical interest. In practice, we find that terminal shocks of pul- 
sar winds have a magnetisation level that is too high to allow for 
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the amplification of short scale electromagnetic fields, so that par- 
ticle acceleration must be inhibited. We have found that gamma- 
ray burst external shock waves propagating into a typical interstel- 
lar medium should lead to strong amplification of the magnetic 
field and to Fermi cycles, even at high Lorentz factor. The ener- 
gies reached by the suprathermal electrons can easily explain the 
afterglow emission through jitter radiation (Medvedev 2000). How- 
ever, if the shock wave propagates in a stellar wind, the upstream 
magnetisation may be too large to allow for particle acceleration, 
eventhough magnetic field amplification should take place. 
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